Streamer Propagation as a Pattern Formation Problem: Planar Fronts 
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Streamers often constitute the first stage of dielectric breakdown in strong electric fields: a nonlin- 
ear ionization wave transforms a non-ionized medium into a weakly ionized nonequilibrium plasma. 
New understanding of this old phenomenon can be gained through modern concepts of (interfacial) 
pattern formation. As a first step towards an effective interface description, we determine the front 
width, solve the selection problem for planar fronts and calculate their properties. Our results are 
in good agreement with many features of recent three-dimensional numerical simulations. 
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Transient discharges occur in various forms |l|, e.g., as 
leaders in spark formation or as streamers in ac silent 
discharges J2J. A common feature is the creation of a 
nonequilibrium plasma through the propagation of a non- 
linear ionization wave into a previously non-ionized re- 
gion. Although it is well known that, depending on the 
polarity of the field, discharge patterns on a larger scale 
may either be fractal [p| , or form more regular non- fractal 
patterns 0, ionization fronts do not seem to have been 
analyzed before as a pattern forming system on scales 
resolving their internal structure. While the idea of a 
shock front or a thin ionization sheet has been formu- 
lated in the literature on streamers in the seventies ||, 
the analytical treatment then frequently was based on 
ad- hoc assumptions and on equilibrium concepts, e.g., 
on the assumption that the high electric field would rise 
the electron temperature and that subsequent ionization 
would be thermal. In the last ten years, models in- 
corporating nonequilibrium impact ionization of neutral 
molecules by free electrons have been investigated both 
numerically ||f7| and analytically J8|. Fig. 1(a) shows 
a snapshot from a numerical study by Vitello et al. J7|] 
of the streamer equations, Eqs. (|l])-(Q) below. Here, the 
evolution of the electron and ion densities between two 
planar electrodes with distance 0.5 cm and voltage dif- 
ference 25 kV is integrated forward in time for parameter 
values describing N2 under normal conditions. At time 
t = 0, the electron density was taken nonzero only in a 
small localized region near the upper negative electrode. 
The figure shows the electron density 5.5 ns later. Each 
contour line indicates the increase of the electron density 
by a decade. The lines enclose a finger-like region (the 
body of the streamer), consisting of a non-equilibrium 
plasma; this region rapidly expands downwards towards 
the anode. In the region outside, the gas is essentially 
non-ionized. The fact that the contour lines in the fig- 
ure are very closely and about equidistantly spaced, il- 
lustrates that the electron density within a zone of the 
order of a few /an grows about exponentially by a factor 
of about 10 10 . Since the total charge density is negligible 



before as well as behind the front, streamer dynamics can 
be viewed as the propagation of a thin charged ionization 
sheet separating a non-ionized high field region from an 
ionized electrically screened region. 

Obvious and up to now unanswered questions are: 
What determines the scale of the pattern (e.g., the lat- 
eral width of the finger- like region), its velocity, the field- 
enhancement near the tip, and what effect do the bound- 
ary and initial conditions have? Triggered by the obser- 
vation of interface-like profiles in the simulations flqjTj and 
by the fact that these are precisely the questions that are 
studied in the field of interfacial pattern formation for, 
e.g., dendrites and viscous fingers 0, we show here that 
pattern formation concepts provide a systematic route to 
unravel precisely these aspects: 

(i) For planar fronts, we trace the "great defect" , "the 
inability of the theory to determine a value for the wave 
speed" [5(c)], to the fact that streamers are an exam- 
ple of front propagation into unstable states in virtually 
all models analyzed |||8|. For such problems, it is well- 
known that the velocity cannot be obtained just by ana- 
lyzing uniformly translating fronts using standard meth- 
ods. In the field of pattern formation, the mechanism of 
dynamical front selection has been understood in the last 
decade [|^,[Ll|. In this paper, we show that this allows us 
to derive all essential properties of planar fronts for the 
model of the recent simulations ||f7| . 

(ii) Clearly, an analysis of planar streamer fronts does 
not suffice to explain the global questions of pattern for- 
mation posed above, such as the field enhancement in 
front of the streamer head or the radius of curvature of 
the tip. However, both the simulations IfljTl] and our anal- 
ysis show that the propagating charge sheet is only a few 
fj,ra thick, while the tip radius and the electrode spacing 
are of order mm or more. This separation of scales, that 
makes simulations so demanding, can be made into an 
analytical tool. Much of our present knowledge about 
similar problems like combustion fronts [ |L2L thermal 
plumes pi and chemical waves H], etc., is based on 
an effective interface description. Such a physically ap- 



pealing formulation can be systematically derived in a 
matched asymptotic expansion to lowest order in the ra- 
tio tin I lout: where the inner length scale 4„ is the thick- 
ness of the front (here the thickness of the charge sheet), 
and l ou t the scale of the pattern, e.g., the tip radius. 
In the effective interface approach that we propose for 
streamers, the charge sheet can be viewed as a weakly 
curved locally almost planar front, since the thickness of 
the charge sheet is much smaller than its radius of cur- 
vature. Like in the other problems, the importance of 
our planar front analysis therefore lies in the fact that, 
apart from curvature corrections, it provides a complete 
solution of the so-called inner problem. 

(Hi) In the non-ionized region outside the streamer, 
the electrical potential <& obeys the Laplace equation, 
V 2 <& = 0. Moreover, our analysis shows that the normal 
velocity of a negatively charged planar streamer front 
(v* below) is a weakly nonlinear function of the field 
E + = — V$ just ahead of it. Both features are reminis- 
cent of the equations for other interfacial pattern forming 
problems like dendrites — e.g., the enhanced diffusion in 
front of a dendrite tip is analogous to the field enhance- 
ment in front of a streamer. Streamers will therefore be 
amenable to the same type of analysis p, p^^3[ . Physi- 
cally, we expect that the interface equations will take the 
form of a conservation equation for a charge sheet (in- 
volving transport terms along the sheet, a stretch term 
due to interface curvature and a term associated with 
charge transport from the plasma behind), supplemented 
with an equation for the front speed that includes curva- 
ture corrections, and an equation for the degree of ion- 
ization created by the front which is not determined by 
any conservation law. The derivation of the appropriate 
equations is left to the future, as the analysis is far from 
trivial due to the coupling to the dynamics of the plasma, 
the fact that the electric field is typically not normal to 
the front, and the fact, that in this fully nonequilibrium 
situation, the curvature corrections do not follow from 
simple thermodynamic considerations. 

We now sketch our analysis ||l5| of planar fronts in 
the streamer model equations fflq) that also underly 
Fig. 1(a). The electron and ion densities n e , n+, and 
the electric field £ obey the balance equations 

dt n e + V R • j e = \n e[ i e £\ a e - E °^ , (1) 

dt n+ + Vr • j+ = \n ef i e £\ a Q e' E °^ , (2) 



and the Poisson equation 



(3) 



V R • £ = — (n+ - n e ) . 
The electron and ion current densities j e and j+ are 

je = -n e He £ - De Vr n e , j+ = , (4) 

so that j e is the sum of a drift and a diffusion term, while 
the ion current j + is neglected, since the ions are much 



less mobile than the electrons. The r.h.s. of Eqs. (1) and 
(2) is a source term due to the ionization reaction: In 
high fields free electrons can generate free electrons and 
ions by impact on neutral molecules. The source term is 
given by the magnitude of the electron drift current times 
the target density times the effective ionization cross sec- 
tion; the rate constant ao has the dimension of an inverse 
length. The exponential function expresses, that only in 
high fields electrons have a nonnegligible probability to 
collect the ionization energy between collisions. 

To identify the proper parameters for the behavior on 
the inner front scale, we note that in the simulations the 
fields just ahead of the front are of order of the threshold 
field E = 2 ■ 10 5 V/cm in Eqs. (Q)-(|). The larger the 
rate parameter ao , the more rapid the impact ionization 
will be, and the thinner the front region. The natural 
length scale for the width of the front will indeed turn 
out to be a^ , which is about 2.3 ^m in the simulations 
|g,0. As the drift velocity of electrons in a field of or- 
der Eq is /i e -Eo, the natural timescale for the motion of 
fronts is then t = (aoMe^o) -1 (~ 3 • 10~ 12 s in ^,Q) and 
the natural scale for the charge density is go = So^o^o 
(« 4.7 10 14 e/cm 3 in [§§). 

For analyzing planar fronts, we now introduce dimen- 
sionless variables, x = Xchq, t = t/to, E = £/Eq, the 
electron density a — en e /qo, and the total charge den- 
sity q = (n+ — n e )e/qo. In these units, the only remain- 
ing dimensionless parameter is the dimensionless diffu- 
sion coefficient D = D e ao/ fj, e Eo. In the simulations 
for N 2 |§,0, this value is about 0.1; for typical gases, 
D ranges from 0.1 to 0.3 Jig] , In these variables, the 
charge conservation equation becomes from (|]), (0), (Q): 
d T q + d x (aE + Dd x <r) = 0. Upon combining this with the 
Poisson equation d x E = q and integrating, we obtain 



d T E 



a E 



D8 X 



(5) 



Here the integration constant is zero because on the inner 
time and spatial scale the charge and electron densities 
vanish for x — > oo, while E(x — > oc) = E + time indepen- 
dent. Eq. (H) and the equation for the electron density 



d T o- = d x {<7 E) + D d 2 x <r + a \E\ e 



-1/\E\ 



(6) 



together constitute the one-dimensional streamer equa- 
tions. These equations have two important classes of 
steady state solutions: the ones with a = 0, E + arbi- 
trary, correspond to the non-ionized state of the gas into 
which the front propagates. The ones with a constant 
(denoted cr~) and E = correspond to the screened 
ionized state behind the front. It is straightforward to 
analyze the linear stability of these states with Fourier 
modes of the form e UJT+lkx . Physically, one expects the 
non-ionized (+) state to be unstable: Any small electron 
density drifts in the field E + and gets amplified due to 
impact ionization while the stabilization due to diffusion 



dominates only at short wavelengths. The corresponding 
dispersion relation uj + = ikE + + \E + \ e~ l /\ E I — Dk 2 
confirms the long wave length instability. It is easily 
checked that screening stabilizes the ionized (— ) states 
at all wavelengths, and that lo~ = — o~ — Dk 2 . 

Propagating streamer fronts are therefore an example 
of front propagation into an unstable state. We thus fol- 
low the common path for such problems 1 1 III]: 

(a) As usual, one can demonstrate the existence of a con- 
tinuous family of uniformly translating front solutions of 
the form o~(£,) and E(£) with £ = x — vt, parametrized 
by the velocity v. This is done by formulating the equa- 
tions for er(£) and E(£) as a flow in the phase space 
(a,E,a') with £ playing the role of a time-like vari- 
able. A front profile then corresponds to a trajectory 
connecting one (— ) = (cr~,0, 0) fixed point with one 
(+) = (0, E + 7 0) fixed point, and the existence and mul- 
tiplicity of these can be studied with counting arguments 
[pLl[| . The family of solutions can be obtained explicitly 
for D = by writing Eq. Q5|) as vd^ In \E\ — er, by insert- 
ing this form into Eq. (|6|) and integrating: we then get 

a[E] = v/(v + E) fL dx exp(— 1/x). This determines 
the flow in phase space for D = 0. 

(b) Physically acceptable front solutions must satisfy the 
additional constraint that the number densities n e of elec- 
trons and n + of ions be positive, i.e., cr(£) > for all £. 

(c) We can show that the condition (b) entails a lower 
bound on the range of velocities. More precisely, one can 
show |ljj that the velocity of physically admissible front 
solutions obey 



V > Vf 



max \v , v 



f ]>0, 



(7) 



where v^ is the fastest nonlinear front pd| ] if it exists. 
Nonlinear fronts correspond to strongly heteroclinic or- 
bits in phase space: they reach the (+) fixed point along 
the eigendirection with the fastest contraction. The ve- 
locity v* = -E + + 2^/D \E+\ exp(-l/|£+|) is the value 
of the velocity v below which the eigenvalues describing 
the flow close to the (+) fixed point become complex, so 
that the cr(£) profiles violate (b) as they oscillate around 
zero far ahead of the front. 

(d) Existing knowledge of front propagation |l(J,|ll| leads 
us to conjecture the following mechanism of front se- 
lection: Fronts emerging from sufficiently localized ini- 
tial conditions ||l6f converge asymptotically to the slowest 
physically acceptable front solution Vf defined in (Q) p7[ . 

We have investigated the existence of nonlinear (tv) 
fronts analytically and numerically and checked the 
above conjecture about dynamical selection by direct 
numerical integration of Eqs. (pf) and (ph. Both qual- 
itatively and quantitatively, our predictions reflect the 
strong asymmetry between fronts moving parallel and 
antiparallel to the field: 

Fronts propagating parallel to the electron drift, i.e., 
into a field E + < 0, are negatively charged. Numerically 



we find no v^ front solutions so that we predict the se- 
lected front velocity to be always the value v* given under 
(c). Here diffusion and ionization help to raise the front 
velocity to a value somewhat larger than the electron 
drift velocity — E + . The degree of ionization a~ behind 
the front only weakly depends on D. The analytic result 
cr~ = a\E = 0] for D = (see formula under (a)) M is 
independent of v and a good approximation for all physi- 
cal values of D, as the lower panel of Fig. 1(6) illustrates. 
Moreover, the values of a~(E + ) extracted from the full 
3D simulations of Vitello et al. [pj (crosses) are close to 
the values we calculate for planar streamer fronts. 

Fronts screening a field E + > are positively charged. 
They can propagate only, if diffusion overcomes the drift. 
As a result, for small D propagating fronts are extremely 
steep and slow. The front velocity vanishes like D, while 
both, the spatial decay rate and the degree of ionization 
behind the front scale like 1/D. In the limit D — > 
this singular behavior can be derived analytically |p.5| . 
For general D, we have predicted the front velocities 
Vf(E + ,D) numerically. They are shown in Fig. 1(b). 

The numerical integration of the initial value problem 
fully supports all our predictions on the asymptotically 
approached front for sufficiently localized initial condi- 
tions. As an example, Fig. 2 shows the spatio-temporal 
development of electron density (a) and field (b) of an ini- 
tial state with E = — 1 and a small charge- neutral Gaus- 
sian ionization seed. The diffusion constant is D = 0.1, 
and the field far from the ionized region is held constant. 
The ionized region initially grows exponentially and the 
electrons drift with the field, till field screening in the 
middle sets in. Then a negative front emerges to the 
right and asymptotically (after At « 20) approaches the 
v* (=1.38) front with o~ = 0.144. The positive front on 
the left initially recedes and then gets stuck by the com- 
bined action of drift and screening. This structure keeps 
slowly evolving in time, however, till after a time of order 
4000, the predicted positive front with ut = 0.0146 and 
cr~ = 6.23 emerges (not shown). 

In summary, we have solved the planar streamer front 
problem. Based on these results, we advocate that 
one should understand streamer dynamics as a two-scale 
problem: on the inner scale, we have a moving ioniza- 
tion sheet, whose thickness « 10/xm is set by the ioniza- 
tion length 1/ao- This interface plays the role of a free 
boundary for the outer dynamics, whose scale is set by 
the global geometry. It is on this scale, that the pattern 
formation problem should be studied. The similarity with 
other well-known interfacial pattern forming problems (a 
Laplace equation for the potential $ in the non-ionized 
region with, apart from curvature corrections, a normal 
front velocity a function of V$), gives us confidence that 
properties of streamer patterns like field enhancement at 
the tip, velocity and tip radius, can be obtained in an 
analogous way P| by including the curvature corrections 
in the resulting effective interface equations. 
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FIG. 1. (a) Electron density profile in a negative 
streamer from the 3D cylindersymmetric numerical sim- 
ulations § of Eqs. (l)-(4). Courtesy of P.A. Vitello. (b) 
Our predictions for planar fronts. Upper panel: v'/D 
(solid) and v* / D (dashed lines) as a function of D for 
positive fronts and for E + = 0.4, 0.6, 0.8, 1.0, 1.2, and 
1.4, from bottom to top. Lower panel: Electron den- 
sity o~ — n e e/qo behind a negative front as a function 
of the field E + before the front for D = (solid line), 
1 (dashes), 3 (dots). Crosses: values of a~(E + ) on the 
symmetry axis in the 3D simulations |7J at times 4.75 ns 
and 5.5 ns, with E + the value of the outer field extrap- 
olated towards the tip, in accord with the asymptotic 
matching prescription. 



FIG. 2. Numerical integration of Eqs. (5) - (6) for 
D = 0.1 in a constant background field E — — 1. Initial 
state at t — 0: lowest line. Each new line corresponds to 
a time step At = 5 and the upper line to t — 100. (a) 
electron density, (b) electric field. 
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